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An antenna array devoted to the autonomous radio-detection of high energy 



cn ■ cosmic rays is being deployed on the site of the 21 cm array radio telescope 

. . in Xinjiang, China. Thanks in particular to the very good electromagnetic 

L^ ! environment of this remote experimental site, self-triggering on extensive 

O ! Ei'ir showers induced by cosmic rays has been achieved with a small scale 

prototype of the foreseen antenna array. We give here a detailed description 

of the detector and present the first detection of extensive air showers with 

^ ■ this prototype. 
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1. Introduction 

The Tianshan Radio Experiment for Neutrino Detection (TREND) is a 
Sino-French collaboration proposing to build a large radio array in order to 
search for Ultra High Energy (UHE) neutrinos. These neutrinos may un- 
dergo charged current interactions with the matter below the Earth surface. 
In order to detect these particles, it has been proposed |l|, |2[ y, \4, la, la, LZf 
to search for very inclined or Earth-skimming extended air showers (EAS) 
induced by the decay of tau leptons produced by charged current neutrino 
interactions in the Earth. The motivation and status of the radio-detection of 
ultra-high energy cosmic rays (UHECRs) and UHE neutrinos will be briefly 
discussed in [section 21 We will state in [section 31 the very good electromag- 
netic conditions at the TREND site and present the proof of principle that the 
TREND antenna array can detect EAS in a stand-alone mode. This result 
constitutes a first and important milestone for the purpose of UHE neutrino 
detection. Given the results and performances obtained with the TREND 
radio-detector, the motivations for a UHE neutrino search with TREND are 
evoked in IsectioiT^ 

2. Development of the radio technique in the context of UHECR 
and neutrino detection 

2.1. Radio antennas for the detection of UHECRs 

Despite intense experimental effort over the last twenty years, the nature 
and origin of cosmic rays above the knee (~ 3 x 10^^ eV) are still uncertain. 
One of the major experimental challenges is the very low flux of high energy 
cosmic rays, thus requiring the use of very large detection surfaces or vol- 
umes in order to acquire statistically significant results |8|. Very encouraging 
results recently obtained by the Pierre Auger Observatory |9|, |lO| indicate the 



relevance of hybrid detection systems deployed over extremely large surfaces 
for the study of UHECRs spectrum, composition and origin. 

Parallel to this, recent works by the LOPES and CODALEMA exper- 
iments have produced interesting results for the characterization of cosmic 
rays above 10^^ eV, using radio antenna arrays triggered by ground detectors 



of EAS [Ul, Il2|, |13|, uM, ll5|, Il6[. These results indicate that radio detection 
could become an alternative or complementary technique to the systems used 
presently for the detection of UHECRs. One of the attractive aspects of ra- 
dio detection is the very low cost of antennas and its easiness of deployment 



over large areas. However, the potential of the radio-detection technique can 
be fully exploited only if working as a self-triggering system. We present in 
this paper the proof of principle for such a self-triggering mode for a radio 
detector of EAS, validated on a small-scale prototype of the TREND setup. 

2.2. Hunting UHE neutrinos with radio antennas 

Radio emission by particle showers was first predicted by Askaryan in 



the early 60's [17|. The use of the radio technique for the search of neutrino- 
induced showers in a natural dense medium was proposed [18|, ll9j in the early 
80's, following Askaryan's suggestions. Latter in the 90's, several projects set 
up antenna arrays to search for radio waves from neutrino-induced showers in 
natural dense medium such as ice, rock salt formations or lunar regolith (see 



for instance |20l . l2ll . l22|, |23|). Even if some very exciting results have recently 



been released (see for instance 2J]), only limits on the UHE neutrino fluxes 
have been set so far by these experiments. Though the attenuation length 
of radio waves in solids are large enough (typically a few hundred meters) to 
equip large detection volumes, the achievement of a better sensitivity requires 
an increase in detection surface and number of (deep)-embedded detectors, 
which could be expensive and time consuming, in the case for example of 
Antarctic ice experiments. 

The setup of radio antenna arrays at the surface of the Earth looking for 
very inclined EAS may represent for UHE neutrino detection a cost-effective 
alternative both to these projects and to the well-advanced EAS ground and 



fluorescence detection techniques also used for this search [25|, |26| . Thanks 
to the important breakthroughs ll|, |l2| performed in air-shower radio detec- 
tion, the idea of setting antennas detecting radio emission in the atmosphere 



has triggered intense interest lately [27|, |28|, 129| . Yet, no complete evaluation 
of this technique's potential has been performed so far. One of the main 
challenge in identifying nearly horizontal showers lies in the existence of an- 
thropic noise sources close to the horizon, together with a possible influence 
of distant storms on the atmosphere electromagnetic condition. A detailed 
analysis of the sensitivity of the TREND antennas to nearly horizontal show- 
ers is beyond the scope of this article. Nevertheless, we will state in this paper 
the uniqueness of the electromagnetic environment of the TREND site, re- 
motely located and surrounded by very close and high (up to 5000 m a.s.l.) 
mountains. 



3. The TREND detection setup and analysis procedure 

The TREND project uses a large part of the existing infrastructures of 
the 21 cm array (21CMA) detector. 

3.1. The 21 CM A interferometer 

The 21 cm array is a radio-interferometer aiming at the study of the 



epoch of reionization [30[. It was completed in 2007 by the National Astro- 
nomical Observatories of China in the Tianshan mountain range (Xinjiang, 
China) |3l| , at an altitude of 2650 m. 

Due to its remote location, the 21CMA site benefits from a very clean 
electromagnetic environment, nearly free of stable sources with significant 
emission beyond the galactic background noise. Only two close-by 10 kV 
power lines and a nearby railway are localized background sources. The 



quality of the electromagnetic environment is illustrated in Figure 1 , which 
represents the radio background observed with a 21CMA antenna directly 
plugged into a 48 dB low-noise amplification system. This spectrum was 
recorded at the galactic signal maximum (18h00 local sideral time). For the 
purpose of illustration, the simulated response of this system to the galactic 



emission is also plotted in Figure 1 



The 21CMA is composed of 10160 antennas, grouped in 80 pods of 127 
antennas each. Commercial log-periodical antennas have been used. They 



consist of 18 side-by-side parallel dipoles (see [Figure 2| ) with geometrical 
patterns designed to work in the 50-200 MHz bandwidth -which corresponds 
to the range of the 21 cm hydrogen emission line for the epoch of reioniza- 



tion 



30[- and achieving antenna gains of ~10 dBi. The dipoles are positioned 
horizontally, resulting in a polarization along the East- West direction, as the 
antennas are orientated towards North. An inclination of 47° with respect 
to ground (so that antennas point towards the Polar Star) corresponds to a 
maximum of the gain pattern for ~40-60° zenith angle over the 50-200 MHz 
frequency band. Over this frequency range, the zenithal beamwidth at 3 dB 
varies between 20 and 60°, while the other antenna characteristics are those 
of standard log-periodical antennas. These antennas do not have the required 
sensitivity towards the horizon for the foreseen search for nearly horizontal 
showers induced by r leptons (see [section 4p . They were nevertheless success- 
fully used for the present stage of the TREND experiment, which consists in 
the validation of the autonomous triggering mode of the radio array on EAS 
induced by cosmic rays. 
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Figure 1; Noise power spectrum measured with a 21 CM A antenna and a 48 dB low noise 
amplification (in black). The green curve shows the spectrum recorded when the input 
cable is disconnected from the antenna. The dashed blue curve is the simulated response 
of an ideal system to the galactic signal only. The frequency range of this measurement is 
0-100 MHz. Few contributions beyond the galactic emission are found above 25 MHz. 



The 21CMA pods form two perpendicular baselines orientated along the 
North-South and East- West directions. They lie in the bottom of two high- 
altitude valleys and extend for 4.0 and 2.8 km respectively (see Figure 3). 



For each pod, the sum of the analog signals from the 127 antennas is contin- 
uously fed into an optical transmitter placed in the middle of the pod. Each 
transmitter is connected through an optical fiber to an acquisition room sit- 
uated in the center of the East- West arm. The signals from the 80 fibers are 
then digitized in parallel by 8-bits ADCs at a sampling rate of 200 MSam- 
ples/s. The data of each pod are finally buffered on 200 MB disks, where the 
signal processing is performed. 

3.2. The TREND prototype setup 

The driving concept of the TREND experiment is that the 21CMA can 
be used as a self-triggered detector for EAS through minor adaptations of the 
existing setup. These adaptations mainly consist in using a single antenna 
instead of a group of 127 as the unit detector, and setting up an on-line soft- 
ware program triggering the writing of the data to the disk independently 
for each channel. 




Figure 2: Left : picture of a 21CMA log-periodical antenna used in the TREND prototype 
setup. Right: picture of several log-periodical antennas and pods from the South-North arm. 
The size of a pod is 20 mx 32 m. 
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Figure 3: The 21 CM A detector layout. Pods positions are indicated with white diamonds. 
The color code gives the altitude a.s.l. in kilometers. The x and y axis correspond respec- 
tively to the West-East and South-North directions. The cross-point of the two baselines 
is taken as the origin of the referential. 



After preliminary tests carried out in 2008, a prototype composed of six 
21CMA-type log-periodical antennas was deployed in January 2009 in order 
to validate this concept. It was tested for several months with various an- 
tenna layouts. The antennas were oriented towards North, thus polarized 
along the East- West axis. Each antenna signal is fed into an optical trans- 
mitter placed in the middle of the associated pod after a 64 dB low-noise 
amplification and a 50-100 MHz filtering. As for 21CMA, the signals are 
sampled by 200 MSamples/s ADCs running in parallel in the acquisition 
room after being transferred through optical fibers. It should be pointed 
out that, as two optical fibers are connecting a pod to the acquisition room, 
TREND can run in parallel with 21CMA without any interference between 
the two experiments. 

Following the transient identification method presented in 'l2i], a software 
trigger is set for each digitized signal, and a waveform composed of 2048 sam- 
ples (10.24 fis) is written to disk whenever a sample value exceeds a threshold 
set as a multiple A^ of a, the instantaneous noise level of the antenna, a is 
calculated every second over 2048 consecutive samples of the buffered data. 
The multiplicative coefficient A^ is chosen at the beginning of the data ac- 
quisition such that the antenna triggers at a rate around 1 Hz. A typical 
value for A^ is 6.5, for which the expected trigger rate is equal to 0.2 Hz for 
a sampling frequency of 200 MSamples/s, assuming a stationary Gaussian 
noise distribution. During quiet periods, we indeed measured trigger rates 
below 1 Hz for A^ = 6.5. 

3.3. Antenna sensitivity 

A 5 ms subset of data is recorded every 5 minutes on all antennas in 
order to monitor the system. The evolution of the noise level of the antennas 
in particular gives a good diagnostic of their sensitivity. The CODALEMA 
experiment has indeed shown that the noise level on a radio antenna should 
follow a periodical evolution, with a higher value when the galactic plane is 
visible in the sky, as the Milky Way is the dominant emitter in the sky at 



radio frequencies [32 



This feature is clearly observed on all TREND prototype antennas, with- 



out any correction nor time averaging (see [Figure 4] ), proving the outstand- 
ing quality of the 21CMA radio environment. Moreover, having identified 
a common physical source shining uniformly over the whole array allows in 
principle to calibrate the antennas responses. This method has not been 
used so far, and the amplitudes of the pulses are presently normalized by 
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Figure 4: Left: evolution of the noise level (in units of least significant bit) for one 
TREND antenna over a period of 14 days. Each point corresponds to the signal standard 
deviation measured over a period of 5 ms in the frequency range 50-100 MHz. Right: 
same data, but with time expressed as local sideral time. The period of high noise level 
corresponds to the crossing of the visible sky by the Galactic plane. Superimposed is the 
expected signal for an antenna pointing to the North, simulated as done in fSSi l. The slight 
mismatch in time is explained by a modest tilt of the antenna from full North towards 
West. 

applying to their raw value a calibration factor equal to 1/cr, a being here 
again the instantaneous antenna noise level measured at the end of the ac- 
quisition chain. As the electronic chain is not the dominant noise source, 
0" is a reliable measure of the antenna sensitivity to electromagnetic signals. 
This method provides a quick way of performing a relative calibration of the 
signals amplitudes. 

3.4- Selection of coincident antenna triggers 

The first step of the off-line data analysis consists in the search for coin- 
cident triggers between several antennas. This is done by ordering in time 
all triggers from a run, and selecting in this list all consecutive values which 
may be associated with the same electromagnetic source. Trigger times for 
antennas i and j are causally linked if they follow the condition: 



\ti 



t,\<^xT, 
c 



(1) 



dij being the distance between these two antennas and c the velocity of light. 
Trigger times ti and tj are defined here as the times when the 2048-sample 
waveforms reach their absolute maximum on antennas i and j respectively. 
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Figure 5: Left: one of the ground layouts used for the 6- antennas TREND prototype 
represented in the cartesian referential defined in \Figure j[ The positions of the antennas 
are indicated by blue triangles, and are labeled by their pod number. The positions of the 
pods are indicated by hexagons. Also plotted is the reconstructed position of the source 
for I642 events recorded during a static source calibration run fsee \subsection 37S\) . They 
are shown as blue crosses, black squares and red circles for events triggering 4-j 5 and 
6 antennas respectively. The true source position is indicated with a yellow star. Right: 
reconstructed positions of the source for the same data, after the cross-correlation treatment 
described in \subsection 3.4\ has been applied. 



The times are corrected for delays associated with the signal propagation 
through the optical fibers and the cables. T is a factor introduced to al- 
low for possible discrepancies between tj and the 'true' trigger time (that is, 
the actual time at which the electromagnetic wave touches antenna i). A 
value T=l.l is chosen in this analysis, which is a safe factor considering the 
TREND timing resolution (see [subsection 3.6P . 



If at least four consecutive triggers from four different antennas follow 



the criterion defined in Equation 1, we consider that these triggers are in 
coincidence and form an event of multiplicity L. We then apply to the 
corresponding waveforms a cross-correlation treatment. Its principle is to 
determine, for the L(L-l)/2 pairs of signals Si and Sj present in the event, 
the delay Tjj which maximizes the cross- correlation coefficient Fy-: 



J- ijVi 



Vj 



IS: 



'xJS,' 



Si{t)Sj(t-Tij)dt, 



(2) 



where Si and Sj actually correspond to a subset of 200 samples around ti 



and tj respectively. For signals associated with close-by sources (for example 
a train moving on the railroad, or a car standing in the middle of the array 



with its engine running as in Figure 5), it appeared that this treatment was 



giving the best results when the signal envelopes were used rather than the 
raw waveforms. In practice however, treatment speed optimization led us to 
perform the cross-correlation treatment on the absolute values of the signals 
rather than the envelopes, the performances of the former being only slightly 
worse. By solving the system of L(L-l)/2 equations t*—t* = Tij according to 
the least squares method, we finally obtain all corrected values t* for antenna 
trigger times relative to an arbitrary reference t*o=0. 

It was observed that for sources standing at fixed positions or moving 
along well-defined trajectories, this treatment significantly improves the pre- 
cision of the reconstructed source position, compared to using the raw trigger 



times (see Figure 5 ). An explanation could be that the L trigger times are 
deduced from an over-constrained system of L{L-l)/2 equations. It can also 
be pointed out that the complete waveform information is taken into account 
in the cross- correlation treatment, while only one sample (the one with max- 
imum amplitude) is used in the determination of the raw trigger time. As 
noise and signal have comparable amplitudes, it is not surprising that the 
later method is less precise. Note nevertheless that the difference between 
the two treatments is smaller than 3 samples (15 ns) for most signals. Given 
its excellent results, the cross-correlation treatment was included in the data 
analysis chain. It should be kept in mind that the waveforms associated with 
these sources are usually characterized by rather long time extensions (sev- 
eral hundreds of nanoseconds). The effect of the cross-correlation treatment 
may therefore be not as significant in the case of prompt signals such as the 
ones expected for the electromagnetic emissions associated with EAS. 

3.5. Reconstruction algorithms 

The following step in the analysis chain consists in the reconstruction of 
the direction of origin of the events formed by coincident triggers. 

Let us first consider the case of two antennas triggering on a wave propa- 
gating spherically from a point source. Their trigger times difference is noted 
Atsph- If we perform the reconstruction of the signal direction of origin as- 
suming a plane wave hypothesis, the resulting trigger time difference Atpian 
will be wrong by a quantity e = Atpian — Atsph given at leading order in d/ R 



10 



by: 

e=^^cosesin2 5, (3) 

where c is the velocity of hght, R the distance to the source and 9 the angle 
between the source direction and d, the vector joining the two antennas. 
Taking d = 250 m (the maximum extension of the TREND prototype array, 
see Figure 5 ) , we find that e <10 ns (our estimated experimental resolution on 
the trigger time measurement, see lsubsection 3.6l) for R>500 m. According to 
this calculation, the discrimination between spherical and plane wavefronts 
is hardly possible for sources further than 500 m with the TREND setup, 
because of its limited extension and modest timing resolution. Monte-Carlo 
studies, carried out with the complete setup geometry, yield similar results. 
Note however that with larger setups, the distance up to which wavefronts 
curvature radii can be reconstructed will increase. 

The maximum development point of a vertical 10^* eV EAS is roughly 
situated at an altitude of 3200 m asl, 500 m above ground at the TREND 
location. The curvature radius of an EAS radio wavefront can therefore be 
reconstructed only for non-inclined EAS with energies above 10^^ eV. This 
most certainly corresponds to a very small fraction of the showers detectable 
by this setup. In the vast majority of cases, the direction of propagation of an 
EAS can therefore safely be reconstructed under the plane wave hypothesis 
with the TREND prototype. This reconstruction is performed through the 
minimization of the following quantity: 

L-l L 
^ plane 



Y.Y.[{x.-x,).k-{t:-t*)\\ (4) 



=1 j>i 

where L is the total number of antennas composing the event, t* and t* are 
the trigger times for antennas i and j corrected through the cross-correlation 
treatment, and Xj, Xj their respective positions. The vector k can be written 
as, k = n/v where n is a vector of norm unity, orthogonal to the wavefront, 
and V is the wave propagation velocity, fixed at the value of the velocity of 
light c. The angular coordinates (zenith angle 9 and azimuth angle 0) of ft 
are the free parameters of the minimization of Fpiane- 

As for local background sources, the signal propagates isotropically from 
the point of emission. In the context of a self-triggering antenna array, we 



11 



therefore found worth of interest to perform the reconstruction of the direc- 
tion of the signal assuming a spherical wavefront. This provides a useful way 
of identifying background sources close to the array. The source location 
and time of emission of a spherical wave can be estimated by minimizing the 
following function: 

where L is here again the event multiplicity, t* is the corrected trigger time 
for antenna i and t^^^ its expected value, which can be written as : 

1 1 \^. Y II 

ti'^ = to + , (6) 

— * — * 

Xi being the antenna position. The source position Xq and signal emission 
time to are free parameters of the fit. As for the plane wavefront case, the 
wave velocity v is taken to be equal to the velocity of light. 

The plane and spherical reconstructions are compared for each recon- 



structed event (see Figure 6). The DMNGB routine from the PORT li- 
brary 33| is used for the minimization of the Fpiane and F^p^ functions. 

A visual indication of the quality of the direction reconstruction can be 
obtained by plotting the measured trigger times versus the values calculated 



with the reconstruction results (see Figure 6 ) . Points deviating significantly 



from the first bisector are an indication of a bad reconstruction of the signal 
direction of origin. A linear fit of this distribution provides a qualitative 
evaluation of the reconstruction. 

3. 6. Reconstruction performances 

The angular resolution of the direction reconstruction can be estimated 
from transient sources crossing the sky. A standard deviation down to 1.5° 
was achieved on a track observed during 3 minutes in one of the TREND runs 
(see [Figure 7\ . Several tracks of this type, most likely due to airplanes, were 
reconstructed in the data recorded with the prototype array, all of them 
yielding similar values for the angular resolution of the plane wave recon- 
struction. It should be noted here that the signals generated by this type 
of sources are likely to be more extended in time than those associated with 
cosmic ray candidates. It is therefore likely that the trigger time determina- 
tion is not as precise for these anthropic sources. In this respect, the values 
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Figure 6: Delay plots for 2 TREND events. Delay plots are defined as the measured trigger 
times versus the values inferred from plane (left) and spherical (right) reconstructions. The 
parameters 9 and (f> are the zenith and azimuth angles of the signal direction of origin. In 
the case of the spherical reconstruction, the parameter p corresponds to the distance to the 
center of pod II4. Antennas are referred by a number (110, 113 to 111) corresponding 
to their associated pod number. A resolution of 10 ns is assumed for the trigger time 
(see lsubsection 3^) . The first triggering antenna is used as the time reference. A linear 
fit of the points is performed. The event represented in the two upper plots was recorded 
during the calibration run performed with a static source placed in the middle of the array 
(yellow star in \Figure 5j) . The source position reconstructed with the spherical treatment 
is [2320.3 m, -62.5 m, 2667.6 m] in the cartesian referential defined in \Figure S\ while the 
true source position is [2322.1 m -62.0 m 2671.2 m]. Obviously for sources so close to the 
antennas, the plane reconstruction fails. 
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Figure 7: Left: directions of arrival for all events detected within 3 minutes in one of the 
TREND data sets. The radial and angular coordinates are respectively the reconstructed 
zenith and azimuth angles 9 and <j) for the plane wave treatment. Zenith angle values of 
9=20, 40, 60, 80 and lOOP are indicated by labels. The color coding for the events plotted 
is the same as the one introduced in \Figure 5\ These events can obviously be associated 
with a mobile source crossing the sky. Right: histogram of the minimal angular distance to 
the source trajectory (yellow curve on the left plot) for the same events. The distributions 
for 4, 5 and 6-antennas events (top to bottom) are well fitted by Gaussian functions with 
a =7.2, 4-4 0''>T'd 1.5° respectively. 



of the resolution inferred from these tracks can certainly be considered as 
conservative. 

The quality of the spherical wave direction reconstruction was also esti- 
mated from data taken with a source (in practice, a car) placed in the middle 
of the array. The source position was measured with a resolution better than 
2 meters after correcting the trigger times with the cross-correlation treat- 
ment presented in [subsection 3.41 (see Figure 5 ) . Sources outside the array 
were also reconstructed with satisfying precision. The spherical reconstruc- 



tion of the direction of origin of the events presented in Figure 7 for example 
resulted in an angular standard deviation comparable with the one obtained 
with the plane reconstruction. 

In the case of the plane wave treatment, the quality of the direction re- 
construction was compared to theoretical expectations assuming Gaussian 
distributed time errors. The theoretical expectations were computed both 
analytically and by Monte-Carlo. For the analytical error estimate, we as- 
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Figure 8: Theoretical angular resolution for the TREND prototype as a function of the 
zenith angle 9 (left, (j) — 90°j and the azimuth angle (p (right, 9 — 45°J. Black curves 
correspond to all 6 antennas triggering. Red curves represent the resolution achieved for 
events with antennas 113, 114, ^5 and 116 only. A timing resolution Ai = 10 ns is 
assumed for this computation. 



sume leading order error propagation. Details are given in [Appendix A[ For 
the Monte-Carlo estimate, random values were generated for the wave direc- 
tion of propagation. The antennas trigger times were then computed for each 
of these directions, and smeared afterward assuming a Gaussian time reso- 
lution. The plane reconstruction algorithm described in the previous section 
was finally applied to these simulated events, and results were compared to 
the true source parameters. Assuming a timing resolution of At = 10 ns 
(corresponding to 2 samples for the digitized data) both treatments yield 
comparable results for the standard error of the reconstructed direction in 
the plane wavefront hypothesis, ranging from 2° for 6-antennas events up 
to 8° for 4-antennas events, in the worst configuration. These results are in 
good agreement with the resolution measured with real sources. Note how- 
ever that a good angular resolution is achieved for sources associated with 
small zenith angles only, as the zenith angle determination degrades when 



moving towards the horizon (see Figure 8 ). This effect is due to the coplanar 
layout of the antennas and could be of some importance in the perspective 
of neutrino detection, in which case the search focuses on nearly-horizontal 
EAS. 
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3. 7. Background event rejection & cosmic ray candidates search 

Background rejection is obviously one of the major challenges for a self- 
triggering experiment like TREND. Even in remote places such as the 21CMA 
site, the environment is not totally free from electromagnetic radiations : 
electrical transformers, trains or the close-by train station as well as the 
21CMA acquisition room are sources of electromagnetic radiations detected 
by the TREND setup. These events have to be discriminated from the ones 
induced by cosmic rays. Fortunately, both types of events exhibit distinct 
features. 

Most background sources are located at ground level. In the case of the 
TREND site, the identified sources are mostly static or follow well-defined 
trajectories (as it is the case for trains for example). The source intensity 
is often variable, with bursts periods and quiet states otherwise. The signal 
duration is likely to be larger than several hundreds of nanoseconds, and a 
high repetition rate (several pulses within 10 /is) is frequent. Finally, the 
corresponding waves are supposed to propagate spherically, with an ampli- 
tude decreasing as the inverse of the distance from the point source. 

Events associated with cosmic rays, on the contrary, follow a random 
distribution in time and in direction of origin in a first approximation. The 
induced antenna signals are generally expected to be very short (<200 ns) 



transient pulses [12, uM, l3J| , roughly symmetrical with respect to the position 
of the signal maximum, and with a subsequent absence of close neighbors in 
space and time. As discussed already, except for the rare case of vertical 
showers with energies above 10^® eV, the radio wavefront can also be approx- 
imated by a plane with the TREND prototype setup. Finally, experimental 
results indicate that the radio amphtude lateral distribution exhibits on the 



average an exponential decrease J35|, l36 



/li = Aoexp(--i), (7) 

"0 

where Ai is the signal amplitude for antenna z, Aq is the signal amplitude 
along the shower axis, d^ the signal attenuation parameter and di the dis- 
tance between the shower axis and antenna i. Values between one hundred 
and few hundred meters are found experimentally for d^ [uJ, ll3|, ll5|, ll6 . 

Taking into account these various statements, we define a noise rejection 
procedure for the data analysis. It should be stressed that at the present 
stage of the TREND project, the priority is to establish the autonomous 
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Figure 9: Top: two waveforms passing rejection criterion 2. Bottom: 2 waveforms 
rejected with cut 2. Time is expressed here in microseconds. 

detection of cosmic rays. The aim of the selection procedure defined in the 
following is therefore the maximization of the background rejection, even at 
the cost of a reduced cosmic ray detection efficiency. This efficiency will not 
be evaluated here. 

1) We first exclude all periods of time for which more than three events 
are reconstructed within three consecutive minutes. These noisy periods cor- 
respond to high emission states for background sources. The risk of misiden- 
tifying a background event is thus increased during these periods of time, and 
it is safer to restrict the analysis to quiet periods, corresponding to events 
rates lower than this threshold. 

2) As transients associated with EAS are expected to last for 200 ns at 
most, all signals with amplitudes significantly higher than the average an- 
tenna noise level for more than 400 ns are discarded, as well as waveforms 
with several transients in the 10 /is sampled window. This cut is illustrated 



through some examples in Figure 9 



3) The accepted events are required to be well reconstructed. The 
quality of the direction of origin reconstruction for a given event can be 
quantified by the x^/iidf value of the linear fit of the delay plot (expected 
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Figure 10: Top : histogram of the X^/ndf values obtained from the linear fit of the delay 
plots (see \Figure 6p for the plane wave reconstruction of 10000 simulated plane wavefront 
events (squares) and 10000 simulated spherical wavefront events (crosses). In the case 
of spherical waves, the wave source position is randomly chosen on the ground within a 
distance of 500 m from the center of the antenna array. Bottom : y^ /ndf distribution for 
the spherical wave reconstruction of the same events. 



trigger times versus measured values, see [Figure 6p . This was studied with 
Monte-Carlo simulations of events associated with plane and spherical wave- 
fronts. The simulation of events associated with plane wavefronts is described 
in [subsection 3.61 The simulation of events associated with spherical wave- 
fronts proceeds as follows : a source position is randomly chosen within 500 m 
from the center of the TREND prototype. The antennas trigger times are 
then computed assuming a spherical propagation of the signal from the point 
source and -as for the plane wavefront case- the computed trigger times are 
then randomly smeared, assuming a Gaussian distribution with a = 10 ns 
for the timing resolution. 

10000 plane wavefront events and 10000 spherical wavefront events were 
simulated as triggering the six antennas of the prototype. The plane and 
spherical wavefronts reconstruction procedures described in [subsection 3.51 
were then applied to these simulated events. A linear fit was performed on 
their delay plots (see Figure 6 ). As expected, the x^/ndf distributions are 
similar for both types of events for the spherical reconstruction, while the 
plane reconstruction for spherical wavefront events is, on the average, worse 



than for plane wavefront events (see Figure 10). It seems reasonable from 



this result to exclude all events associated with a x^/ndf value above 2 for 
the plane reconstruction. However, data tend to exhibit larger x^/ndf val- 
ues than what is observed with simulations. This is most certainly due to 
the systematic uncertainties affecting the measurement of the signal delay in 
the coaxial cables and optical fibers, estimated at ±0.5 sample (i.e. 2.5 ns) 
for each antenna. The cut set to the plane waveform reconstruction of the 
direction of origin in our cosmic ray search procedure is therefore : 

X^/ndf < 5. 

4) Since most background sources are expected to be located at ground 
level, the signal/noise ratio should increase when moving towards the zenith. 
In order to establish the autonomous radio-detection of cosmic rays, we there- 
fore choose to limit our candidates search to events reconstructed with zenith 
angle values below 65° : 

"plane < 00 . 

5) It has been shown in [subsection 3.5] that the time resolution and array 
extension of the TREND prototype do not allow for the measurement of 
the curvature radius of the EAS radio wavefront, except for the very rare 
case of nearly vertical showers with energies above 10^^ eV. Both plane and 
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spherical wavefronts treatments should therefore yield similar results for the 
reconstructed direction of origin of the signal. We therefore request the 
following conditions for the results of the direction of origin reconstruction : 

y {dsph — dplanY + {'t'sph — <Pplan)^ < dcut, 

where dcut = 10, 7 and 4° for 4, 5 and 6-antennas events respectively. These 
values roughly correspond to the quadratic sum of the average error for the 
spherical and plane wave reconstructions. 

These various selections were applied to simulated data in order to es- 
timate the background rejection efficiency. This Monte Carlo simulation 
was produced following the method described in [subsection 3.6t with ground 
source positions randomly chosen within 10000 m from the array center. Over 
10^ simulated events, none passes cuts 3 to 5. This is an interesting result, 
even though it has a limited validity only, as the simulation performed here 
presents very basic characteristics, far from reality. 

6) The amplitude information should in principle provide a powerful lever 
for background rejection : as pointed out already, the lateral profile of the 
signal amplitude should preferentially decrease exponentially for an event 



associated with an EAS (see Equation 7). The combination of this rapidly 



decreasing signal amplitude with a large curvature radius of the wavefront 
could constitute a distinct signature of cosmic ray candidates as demon- 
strated recently |13|. Indeed, background events should either be associated 
with large curvature radii and roughly constant amplitudes in the case of 
distant sources, or spherical wavefronts and rapidly decreasing amplitudes 
for close sources. It is therefore believed that background rejection efficiency 
of an antenna array will increase with its size. The reduced extension of 
the 6-antennas TREND prototype array (~200 m in its largest dimension) 
and the limited precision of the calibration technique presently used (see 
[subsection 3.31) do not allow for a reliable reconstruction of the shower core 



position and lateral amplitude profile. This complete reconstruction of the 
shower characteristics for events passing cuts 1 to 5 will therefore be per- 
formed (see [subsection 3.8^ . but not included in the cosmic ray candidate 
search procedure. The minimization function for the shower geometry recon- 
struction is assuming an exponential decrease of the lateral amplitude and is 
given by : 



shower 



y^(-So - ao X di - Sif, (8) 
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where the sum runs over all triggering antennas, Oq is given by Oq = j^ (lo)xd ' 
do being the attenuation parameter of the exponential fit of the lateral am- 



plitude profile as defined in Equation 7 Sq is the electromagnetic field am 



plitude in decibels along the shower axis, di the distance from antenna i to 
the shower axis and Sj the measured antenna amplitude in decibels. The 
parameters sq, cto, and the shower core position (xq, yo) are free parameters 
of the minimization. The shower direction is given by the vector n previously 
determined through the plane wave direction reconstruction of the event (see 



Equation 4). As before, the function minimization is performed with the 



DMNGB routine from the PORT hbrary [33 



3.8. Results 

In this section, we apply the selection procedure defined in the previous 
section to the data recorded in 2009 with the 6-antennas TREND prototype. 

The setup ran for 584.7 live hours in 2009. The noisy periods amount 
to 31% of this total. In the 403.4 hours selected as quiet periods, 2275 



events are successfully reconstructed (see [Figure 11[ ). The vast majority of 
these events lie along the horizon and are mostly associated with background 
sources. It is noticeable that few of these ground events are reconstructed 
in the North direction. This can be explained by the fact that no potential 
background sources can be found North from the array, while several exist in 
the South (railroad). East (train station) or West (DAQ room and electrical 
transformers) directions. 

After applying cuts 2 to 5 to the 2275 events recorded in quiet periods, 
only 37 events are left. A detailed analysis ends in a further rejection of 12 
of these candidates, which are possibly related to other events reconstructed 
close in time or space. On the contrary, all 25 surviving candidates are 
clearly isolated in time and space, and are therefore considered as cosmic 



ray candidates. We plot in Figure 12 the reconstructed arrival directions for 



these events. It is worth noting that 20 of them lie in the Northern part of 
the sky, while there are significantly less background sources at ground level 
towards North. As various biases have not been taken into account in the 
present analysis (in particular the antenna lobe), it is certainly premature to 
relate the observed North-South asymmetry to a possible geomagnetic origin 
of the radio signal, as done by CODALEMA [lit . 

A complete reconstruction of the shower characteristics was successfully 



performed on 18 of these 25 candidates, as shown in Figure 13 for two nice 



examples. The average value for the attenuation parameter is (rfo) — 200 m, 
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Figure 11; Left : plot of the direction of origins for the 2275 events reconstructed during 
quiet periods. The color coding for the events plotted is the same as the one introduced in 
\Figure 5[ Right: histogram of the reconstructed azimuth (top) and zenith (middle) angles 
for the same events by bins of 1° . Also shown in the bottom plot is the cos 9 distribution 
for these events by bins of 0.02 radians. The direction (j) ^105 ° corresponds to 2 electrical 
transformers close to the 21 CM A acquisition room. 




Figure 12: Repartition of the direction of arrival for the 25 cosmic ray candidates detected 
with the TREND 6-antennas prototype. Zenith angle values of 9=20, 40, 60, 80 and 100^ 
are labeled. The color coding for the events plotted is the same as the one introduced 
in \Figure 5[ Also shown as a green star is the direction of the geomagnetic field at the 
TREND location. 



22 





1113 






a 


I 








x^^ 






u 


\^ 






•a 20 


\ 


J 


- 


3 




X^^ 




1 


110 


114 IV 
115 


n 116 


-." 




^. 


a 






^\^ 


M 






J\ 


GO 






117 1 


10 









3 



73 



I 

a 
cJ5 



-50 50 100 150 200 250 300 

Distance to shower axis [m] 



JU 

45 


\ 115 




- 


40 


V 




- 


35 


\ 


N,"6 


- 


30 






- 


25 


114 


113" 


"^ 


110 


20 




117^ 


v. 


15 









20 40 60 80 100 120 140 160 

Distance to shower axis [m] 



C'ori' position 

\VK = 2251m, SN = 2«in 



e 

Z 



IK)' 



115 117 



114 116 



2100 2200 2300 2400 2500 

W-E [ml 



100 



Z -100 



-200 



-300 



110 



114 117 



113 ^^^ 116 
•115 



Core position 
WE=2322m, SN=-93in 



2200 2300 2400 2500 
W-E [m] 



Figure 13: Left : signal amplitude as a function of antenna distance to the shower 
axis and exponential fit for 2 cosmic ray candidates. Only statistical uncertainties on the 
amplitudes are shown. For the second candidate, the amplitude of antenna 11 ^ is not taken 
into account : as the antenna was placed atop a 5 m high pole at this time, its relative 
calibration is not reliable. Right : projection of the reconstructed shower axis (red line, 
dashed is underground) on the detector plane. The core position is represented by a red 
star. A disk of radius proportional to the signal amplitude is drawn around each antenna 
position. Note that these 2 events were recorded with 2 different antennas layouts. The 
shower plane directions are (9 ^ U ± 8° , (/) = 2 ± 2°) and (0 = 42 ± 5°, (/> = 13 ± 3°; 
respectively. Values do— 4-62 m and do=158 m are found for the exponential fit of the 
lateral amplitude profile (see \Equation T] ). 
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result consistent with calculated values 37 and available radio-detection 



data [I3|,ll6|]. It should be stressed however that in these experiments, ground 



detectors are used to trigger the antenna array. 

3. 9. Validation of the cosmic ray search procedure 

In January 2010, 15 TREND antennas were set up around the cross-point 



of the 21CMA baselines (see Figure 14), covering a total area of 350 mx 800 m 



~ 0.2 km^. Data have been recorded with this setup since then, and the cos- 
mic ray candidates search procedure described in [subsection 3.71 has been 
applied to the reconstructed events. 

An array composed of three scintillator detectors was installed at the same 
location. The three detectors are separated by ~200 m from one another. 
Each of them is composed of a 0.5 mxO.5 mx2 cm standard plastic scintil- 
lator in direct view of a Photonis XP2020 photomultiplier tube (PMT). The 
PMT signal is converted by an optical transmitter working in the 20-200 MHz 
frequency range and sent to the acquisition room through an optical fiber. 
Similarly to an antenna signal, it is then digitized by an 8-bits ADC work- 
ing at a rate of 200 MSamples/s. A triggering procedure similar to the one 
described in [subsection 3.21 is applied to the scintillators signals. The value 
of the factor A^ is set in order to obtain an individual trigger rate around 
25 Hz for each scintillator. A selection of triggers in coincidence on the 3 
scintillators is performed in an off-line analysis procedure. 

The two systems have run in coincidence for 19 live days so far. Dur- 
ing this period, 620 3-fold scintillator events were recorded. The rate of 
random 3-folds coincidences being lower than 3xl0~^/day according to cal- 



culations (see Appendix B ) , these coincident events have to correspond to 
EAS in their vast majority. A reconstruction of the direction of origin was 
performed for these events a ssuming a plane wavefront. The distribution of 
these reconstructed directions is compatible with what is expected for cosmic 



rays (see [Figure 15| ) . 

During this period of hybrid operation, three cosmic ray candidates 
selected in the radio antenna data following the procedure described in 



[subsection 3. 71 were recorded in coincidence with three scintillators (see Figure 16 ) 



and two with two scintillators. We have carefully examined these coincidences 
and the various hypotheses for their origin. 

Firstly, we can exclude the occurrence of random coincidences. The rate 
of random coincidences between two systems A and 5, triggering on uncor- 
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Figure 14: Layout of the TREND setup installed at the cross-point of the two 21 CM A 
baselines. Antennas are shown as blue triangles, scintillators as red squares. Also shown 
as black diamonds are the 21CMA pods. TREND detectors are labeled by numbers corre- 
sponding to the pods they are connected to. 



related stationary noises, is computed in Appendix B It writes as: 



/, 



/a/ 



B 



rdm 



fA + fE 



;i-exp(-(/^ + /B)At)), 



(9) 



with Ja and Jb the trigger frequencies in A and B, and At the time window 
considered for a coincidence. Over the 19 days of data-taking, 620 events 
were detected in coincidence on the 3 scintillators, corresponding to a trig- 
ger rate /a = 3.8 xlO"'* Hz. It has also been observed in practice that Jb = 
10 Hz is a safe upper limit for the rate of coincidences on 4 antennas or more. 



Injecting these values in Equation 9 it is found that the rate of random coin- 
cidences between the 3 scintillators and 4 antennas in a time window of 2 fis 
can be estimated at less than one per year. The chance to observe within 19 
days two random coincidences between two scintillators and four antennas 
or more is also negligible. The five hybrid coincidences detected by TREND 
therefore have to be related to the same physical source. 

Secondly, the influence of PMT radiations on the antennas can also be 
excluded. For four out of these five coincidences indeed, antennas triggered 
before the scintillators. The time ordering of the triggers is also inconsistent 
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Figure 15: Left: repartition of the direction of arrival for the 620 3-fold events recorded 
with the scintillator array. Right: distribution of zenith and azimuth angles for the same 
data by bins of 6 and ^CP respectively. The zenith angle distribution is fitted by a curve 
given by -^ = {a + b9) sin cos x ^/)-eo)/ei , corresponding to the expected zenith distri- 
bution for a scintillator array \l aj . The adjusted values for a, b, Oq and 9i are respectively 
68°~^ , -1, 4-5° and 1S° . These are standard values given the angular resolution of the 
array. 





Radio antennas 


Scintillators 


Event 


e, <p 


e, <p 


A 


52±1, 195±2 


49±3, 191±4 


B 


61±3, 359±2 


67±5, 3±4 


C 


42±1, 55±4 


36±3, 56±5 



Table 1: Reconstructed zenith and azimuth angles (in degrees) for the 3 hybrid events 
detected with the TREND setup (plane wave hypothesis). The direction reconstruction 
has been performed independently for the antenna and ground array data. All results are 
statistically compatible. The uncertainties for the antenna and ground array direction 
reconstructions are estimated using \Equation A.1\ by assuming a time resolution of 10 
and 20 ns for the antennas and scintillators respectively. 
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Figure 16: Waveforms for the hybrid event A, involving 4 antennas and 3 scintillators 
in coincidence. Delays induced by signal propagation through the fiber and cables are not 
taken into account for in this display. 
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with the transit time of an electromagnetic wave from the PMTs to the an- 
tennas in all five cases. 

At last, direction reconstructions were performed separately for the ra- 
dio and scintillator data for the three events with three scintillators. The 
results of these two independent reconstructions are in excellent agreement 
within uncertainties for the three events (see lTable 1|) . This agreement gives 



us strong confidence that these events have indeed been induced by EAS. 
In addition, for the two coincidences with two scintillators, the scintillators 
trigger times are consistent with the values expected from the direction of 
propagation of the wave reconstructed from the radio data. 

These five coincident events therefore constitute an unquestionable proof 
that autonomous radio-detection of cosmic rays has been performed by TREND, 
and that the selection procedure defined in [subsection 3.71 has allowed to sin- 
gle them out from the background events. Apart from some partial measure- 



ments [12| quoted earlier, our work establishes for the first time a reliable 
background rejection method allowing the operation of radio antennas in a 
stand-alone mode. 

4. Outlook 

The ultimate goal of the funded TREND project is the search for signa- 
tures of cosmic neutrinos with energies higher than 10^^-10^^ eV. We think 
that the 21CMA site is very well suited for this search. As pointed out in 
[section 2"| it is indeed surrounded by close and high mountain ranges situated 
~100 km north of a vast depression. This environment results in rock thick- 
nesses of several tens of kilometers and atmosphere lengths of several kilome- 
ters (depending on the observed direction). For energies above 10^^ eV, this 



topology favors the observation of r-induced showers [27|, |28|, |29[. In addi- 



tion to this and its excellent electromagnetic environment (see [subsection 3.11 



and [subsection 3l3|) . the TREND site may also benefit from a very low thun- 



derstorm or lightning activity [39]. These features, together with the low 
cost, easy deployment and maintenance of the TREND detection equipments, 
could be decisive elements for the success of the TREND project. 

Several adaptations are however necessary to complete the TREND setup. 
An antenna design with an optimized sensitivity along the horizon is in par- 
ticular required, as well as the development of appropriate signal identi- 
fication and background rejection procedures. These various elements are 



currently under study and will be presented in a forthcoming paper |38 
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together with an estimation of the TREND sensitivity to UHE neutrinos. 
The purpose of the present paper is to state the excellent electromagnetic 
conditions at the TREND site and demonstrate our ability to perform the 
self-triggering detection of EAS, a prerequisite to the final goal of the Tian- 
shan Radio Experiment for Neutrino Detection project. 

5. Conclusion 

We have given a detailed presentation of the on-going TREND experi- 
ment. We described the analysis procedure used for the data recorded in 
2009 on a 6-antenna prototype working in a stand-alone mode. Using suc- 
cessive rejection criteria to isolate true physical events, we have selected 25 
events which present all characteristics of cosmic-ray induced signals. We 
have shown the ability of our setup and analysis procedure to perform the 
autonomous radio- detect ion of EAS by using a scintillator array running inde- 
pendently in coincidence with an upgraded version of the TREND prototype. 
With this hybrid setup, several cosmic ray candidates selected from the radio 
data were independently detected by the scintillator array. Follo wing the first 



results obtained on small self-triggering radio-antenna set-ups [12|, |40|, |41 
TREND has here proposed and validated an identification method for EAS 
which clearly demonstrates the possible use of a radio-antenna array in a 
stand-alone mode for their investigation. Based on these observations and 
the fulfillment of prerequisites regarding the electromagnetic and topological 
environments of the TREND experiment site, it is now our goal to exploit 
this analysis procedure with an extended antenna array dedicated to the 
identification of UHE neutrinos. 
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Appendix A. Error propagation for the plane wave reconstruction 

Let us assume that we have a set of n detector units located at positions 
Xj, the arrival time of a plane wave at these locations at times tj with i = 
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l,...,n. Let us further assume that the uncertainties on the unit detectors 
positions and on the wave propagation speed can be neglected. We write 
Ati = ti — ti the errors on the arrival times, with t, being the true arrival 
time of the plane wave. Similarly, the angular errors on the plane wave 
direction (parametrized in spherical coordinates) are written A^ = 9 — 9 and 
A0 = — 0. At leading order the error Ak on the wave vector k writes: 

Ak = dek A9 + d^k A(j) + k, (A.l) 

where d^ is the partial derivative operator as x. By definition, ti and k 



cancel each of the terms of the sum in Equation 4 such that, at leading order. 



FpianeiU', 9, 0) = Fpiane{Ati] A9, A0). Furthermore, 9 and are calculated by 
minimizing the plane wave residuals, Fpiane, hence A^ and A0 are related to 
Atj's through: 

deFplane = d^pFpiane = 0, (A. 2) 

which yields a linear set of equations. Inverting the latter equations for 
various tj's one finds the error propagation coefficients, written as: 

dt,.9 = -r{akfl,e^<f>,<f> — Q'kfl,<j>^e,<i)) (^-3) 

dtk4> = -r {(^kfl,4> h,e - akflfi h,^) , (A.4) 



with 





at, 


n— 1 n 


X,) 


■ da;k 


b. 


,y ~ 




0'i,j,x ' 


^ij,y 



(A.5) 
(A.6) 



By extension we write Xq = ^ Xi/n, the middle position of the set of de- 
tector units. The determinant A = be^e b^^^ — b^^^ falls from the inversion of 
the equation set IA.2I 



Appendix A.l. Standard deviation on angular coordinates 

In the particular case of independently distributed time errors of stan- 
dard deviation at,., the resulting standard deviations ag and a^ on 9 and 
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at leading order can be written as a squared sum of the error propagation 
coefficients : 

n n 

^l = Y.^\ef< ' ^l = 11^9t,ct>fal. (A.7) 

k=l k=l 

Let us further define the absolute pointing accuracy as \/{'ip'^), with ip the 
angle between the true wave direction and the reconstructed one. At leading 
order, ip'^ is linearly related to A6'^ and A0^ and the pointing accuracy simply 
goes as: 

{i;')=a',+sm'i9)al. (A.8) 

Appendix A. 2. Statistical significance of the residuals 

Further expanding A9 and A0 at leading order in Atj's, the plane wave 
residuals, Fpiane is: 

Fplane = AT'CAT, (A.9) 

where AT is a size n vector of the Atj's. C is a n x n correlation like matrix, 
whose coefficients are given as: 

Ck,i = dt^edti6hefi + dt^(l)dti(t)h^^^ + {dt^edti(t) + 

dt^^(j)dt,6)be^^-n{akfi,edti0 + ai^o^e dt^^O + (A.IO) 

ak,o,4>dti4> +ai^(},<i,dt^:(j)) +n6k,i-l, 

where 6k^i is the Kronecker 5-function. 

Let us consider the particular case of independent and identically dis- 
tributed {iid) time errors, Atj, following a centered Gaussian law of standard 
deviation at- The matrix C is symmetrical and positive-defined, hence it 
diagonalizes over an ortho normal basis, Ui, yielding : 

n 

with Aj > the eigen values of C. Following from orthonormality, the base 
vectors Ui are also centered Gaussians iid with the same standard devia- 
tion at- Consequently, the normalized plane wave residuals Fpiane/<^t ^^^ 
distributed as linear combinations of x^ variables, with positive coefficients. 
The cumulative density function, Pi Fpinnp < c) can efficiently be computed 
using Ruben's expansion in x^ series [42|, |43| . See for example the algorithm 
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AS 204 4^ for a numeric implementation. This procedure provides the sta- 
tistical significance or p-value P{Fpiane < Fpiane) ^^ ^^^ measured plane wave 
residuals, F"^^^^ when testing the hypothesis {Hq: the wavefront is plane like 
with centered Gaussian iid time errors of standard deviation at}- 

Appendix B. Random coincidences 

Let us first consider the case of a single detection unit triggering on a 
random stationary noise at a rate /. Note that the detection unit might 
be a set of various sub-systems as long as the triggers it provides can be 
considered as stationary. Starting the observation at an arbitrary time to; 
we recall that the elementary probability, dp for the first trigger to occur 
exactly in the time interval [t; t + dt], t > to, goes as: 

dp = fexp{-f{t-to))dt, (B.l) 

while the probability, p, to have no triggers in the time interval [to; t] is given 
as: 

p = exp(-/(t-to)). (B.2) 

Appendix B.l. Random coincidence probability 

Let us now consider the case of n detector units triggering on uncorrelated 
stationary noises at rates fi {i G [!,'",]). Starting the observation at time 
to = let us assume that the following event occurred: {Sq: each receiver 
triggered only once, at time ti , and we have to < ti < t2 < ... < tn}. The 
event So is the sum of n independent events {Si: the i*^ trigger occurs at time 
ti starting from tj_i and there is no other unit triggering in the time interval 



[tj_i;tj]}. Following from Equation B.l and Equation B.2, the elementary 



probability, dpi, for the event Si writes as: 

n 

dpi = fiexp{-fi{ti-ti_i))dti Yl exp(-/j (tj - ti_i)) 

= fi exp(-/s {ti - ti_i)) dti, (B.3) 

where /s = ^"=i fi- Since the event So is the sum of the independent events 
£j's, the elementary probability d"'p for the event So is the product of the 
elementary probabilities dpi, yielding: 

n 

dy = exp{-fj:tn) Hfidt,. (B.4) 
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Integrating the latter probability density over ti,t2, ...,t„_i one gets the ele- 
mentary probability dp for the n^^ trigger to occur at the time tn, as: 

I ...J^ dy=\Jlfij^^f^exp{-f,,Qdt^. (B.5) 



.1=1 



The probability P(t„ < t) that the given succession of triggers has occurred 
at time t is then given by one more integration of dp over [0,t]. It can be 
checked that this probability does not depend on the time ordering of the 
n triggers. Hence, the probability is the same for any permutation of the 
sequence of triggers. It only depends on the time tiast of the last trigger and 
on the rates /». Therefore, considering the n\ permutations of the sequence 
of triggers the probability P(tiast < t) that each of the n units triggered once 
and only once in the time interval At = t — to is n\ P{tn < t), yielding: 

P{tiast <t)=n\(f[pA (1 - T„_i(/s At) exp(-/E At)) . (B.6) 

where pi = fi/ Jt. is the probability for a trigger to occur on the z*^ unit 
and Tn{x) = Yll^=o ^^/^^- is the n^^ order Taylor expansion of the exponential 
function. 

Appendix B.2. Random coincidence rate 

Let us come back to the sequence of n triggers triggering at the times ti < 
^2 < ••• < ^n- Taking the 1** trigger time, ti = 0, as the time origin a similar 



computation as for Equation B.3 and Equation B.4 yields the elementary 
probability d^pi for getting the time sequence t2 < ... < t„, given that no 
other trigger occurred on the 1** unit. It goes as: 

n 

dy, = exp(-/s U) n /^ dU, (B.7) 

i=2 

with /s the sum of all rates, including the rate /i for triggering on the 1** 
receiver. Integrating the latter equation one can compute the probability 
Pi(t„ < t) that this trigger sequence has occurred at the time t. The rate 
Ri(tn < t) of such events is then given by the product of the triggering rate 
/i on the 1'^* unit times the probability Pi(t,„ < t) for the right triggering 
sequence to follow. Once more, one can check that this rate does not depend 
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on the particular ordering of the trigger sequence. Therefore, considering all 
permutations of the trigger sequence, the rate Ritiast < t) of coincidences of 
the n detector units, in a time window At = t — t first, is: 

Ritiast < t) = n\ f^ (f[pA (1 - T„_2(/s At) exp(-/sAt)) , (B.8) 

where it must be understood that it is assumed that each receiver triggered 
once and only once in the time window At. Nevertheless, a similar reasoning 
as used previously holds that the probability to get one more trigger from 
any detector unit is negligible within the time window At provided that the 
product /s At is small enough. 
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